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Abstract. The angular two-point correlation function between background QSOs and 

foreground galaxies induced by gravitational lensing is derived. It is shown that the 

shape of this correlation function depends sensitively on the spectrum of the density 

fluctuations in the Universe, thus providing a possibility to distinguish between different 

models for the spectrum. Using numerical large-scale structure simulations, I estimate 

that the QSO-galaxy correlation function can be measured from galaxy counts down 

^ to ~ 21""* . . . 22"*^ magnitude in fields with radius ^ 25' around 50 . . . 100 QSOs with 

0\ redshift z > 1. Since the QSO-galaxy correlation function is proportional to (—a — 1)6, 

I where b is the biasing factor of galaxy formation and a is the effective slope of the QSO 

,__! number counts, steep number counts are favorable for this kind of analysis. I show that 

^ a = —4 ... — 5 can be achieved with the 1-Jansky sample of radio-loud QSOs when 

QQ the double-waveband magnification bias is employed. Moreover, the cross-correlation 

CN analysis allows to determine the galaxy-formation bias factor h. 
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Q 1 Introduction 

^ It was shown in a sequence of previous papers that there is evidence for highly significant 

correlations on angular scales of several 10 arc minutes between radio-loud, high-redshift 
QSOs and foreground optical galaxies (Fugmann 1992, Bartelmann & Schneider 1993b), 
IRAS galaxies (Bartelmann & Schneider 1994a) and diffuse X-ray emission (Bartelmann, 
Schneider & Hasinger 1994b). In an earlier study, it was demonstrated on the basis of 
numerical large-scale structure models that gravitational lensing by large-scale structures 
provides a viable explanation for such correlations (Bartelmann & Schneider 1993a). The 
basic underlying idea was the following. Dark matter, which is usually assumed to dom- 
inate the mass density in the Universe, must be inhomogeneously distributed in space, 
because otherwise structure formation can hardly be understood. A matter distribution 
with spatially varying density acts as a gravitational lens, magnifying and distorting 
images of background sources. In a flux-limited sample of background sources such as 
high-redshift QSOs, such magnified sources will preferentially be included. This effect. 
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termed magnification bias, is extensively described in the literature (see, amongst nu- 
merous others, Schneider, Ehlers, & Falco 1992, Chap. 12, Narayan & Wallington 1993, 
and references therein). This implies that sources from a flux- limited sample will prefer- 
entially be found in the background of matter overdensities. From the biasing hypothesis 
of galaxy formation (e.g.. Kaiser 1984, Dekel & Rees 1987), it is expected that galaxies 
form preferentially where the dark matter is overdense. Combined with the magniflca- 
tion bias caused by these overdensities to a flux-limited sample, this implies that such 
background sources are preferentially observed where the foreground galaxy density is 
enhanced, and this yields gravitational-lensing induced correlations between background 
QSOs and foreground galaxies. 

The occurrence of such correlations has been tested and established on a signiflcance 
level of up to 99.8% in the previously cited papers. There, Spearman's rank-order corre- 
lation test was employed. This test has the advantage of being robust, parameter-free, 
and highly sensitive (e.g., Kendall & Stuart 1973), but it has the disadvantage that it 
yields a 'binary' result, since it merely answers the question of whether QSOs and galax- 
ies can be considered correlated or not. Nevertheless, the correlation tests performed 
unambiguously demonstrate that such a correlation exists, although they do not allow 
to quantify the strength of the correlation; note that a significant correlation does not 
necessarily imply a strong correlation if a highly sensitive correlation test is used. 

From the arguments stated above, it is clear that, if gravitational lensing is indeed 
the reason for the detected correlations, these prove the biasing hypothesis of galaxy 
formation. If this hypothesis were not at least qualitatively correct, galaxies were not 
correlated with dark-matter overdensities, and then the distribution of dark lenses would 
be random compared to the distribution of galaxies. 

This paper addresses the question of how the biasing hypothesis of galaxy forma- 
tion might be quantifled on the basis of the gravitational-lensing interpretation of the 
correlations between background sources and foreground galaxies. Such a quantiflcation 
is interesting in various respects. First, it constitutes a further observational test for 
such cosmogonic scenarios which are fairly specific in the amount of biasing they require 
to describe the observed galaxy distribution from the theoretically modeled dark-matter 
distribution (e.g.. White et al. 1987, Cen & Ostriker 1992, Ostriker 1993). Second, it 
provides a lower limit for the mass-to-light ratio on cosmological scales beyond the clus- 
ter scale. This is because the bias factor quantifies the ratio of the fiuctuations of the 
galaxy number density relative to the density fiuctuations. Homogeneously distributed 
dark matter does not contribute to the magnification bias due to gravitational lensing 
because it does not magnify background sources, and thus only deviations of the cosmic 
density from its mean contribute to the effect. 

In Sect. 2 of the present paper, I derive the angular cross-correlation function ^qg(0) 
between background sources and foreground galaxies as expected from gravitational lens- 
ing by large-scale structures. This is done in close analogy to Kaiser (1992). I show that 
^qg(^) is proportional to (— a — 1), where a is the effective double-logarithmic slope of the 
luminosity function of the background sources, and to the effective biasing factor b of the 
foreground galaxies on the scale considered. The shape of ^qg(0) depends strongly on 
the shape of the dark-matter perturbation spectrum Ps{k). A series expansion of ^qg(</') 
for <C 1 shows how parameters determining the shape of Ps{k) can be derived from 
^qg(^)- In Sect. 3, I compare the analytic results obtained in the previous section with 
a correlation function derived from numerical simulations incorporating a simple recipe 
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for galaxy formation. It is argued there that in reahty the biasing parameter may be a 
function of scale, as it is in the numerical simulation. This scale-dependence of b slightly 
steepens the correlation function. The numerical simulation is intended to demonstrate 
the feasibility of determining ^qg(</') from real data if sufficiently deep (^ 21*"^ magni- 
tude) galaxy observations in sufficiently large (^ 25') fields around 50 . . . 100 background 
QSOs are available. Bootstrap error estimates are given. In Sect. 4 I argue that the 
multiple- waveband magnification bias first discussed by Borgeest et al. (1991) effectively 
steepens the luminosity function to cumulative slopes of a = —4 ... — 5 in the 1-Jansky 
sample, dependent on the optical flux limit imposed. A slope as steep as that was re- 
quired in earlier theoretical work (Bartelmann & Schneider 1992), and it is shown that 
it can indeed be achieved with real QSO samples. Sect. 4 summarizes the results. 



2 The angular QSO-galaxy correlation function from weak 
lensing 

2.1 Light propagation 

It was shown by Kaiser (1992) that the (two-dimensional) deflection angle of a light 
ray propagating through an Einstein-de Sitter background universe which is slightly 
perturbed by density inhomogeneities with Newtonian gravitational potential <? is given 

by 

2 Z"™ 

aj{e,w) = / dw'{w-w')^j{w'ei,w'e2,w';l-w') . (2.1) 

^ Jo 
To render this paper self-contained, I re-derive this equation in Appendix A. The index on 
^, preceded by a comma, denotes partial differentiation with respect to the coordinate 

— * 

Wj. w is the comoving distance in units of twice the Hubble length, 9 is the angular 
position vector of the light ray on the observer's sky, and the argument 1 — iw of ^ is the 
conformal cosmological time at which the potential perturbation <? is considered. The 
important approximation which entered the derivation of Eq.(2.1) was that the potential 
gradient is integrated along the unperturbed path of the light ray. 

The deformation of an infinitesimal light bundle is determined by the shear tensor 

ljk{9,w) = ^^Jll^ = -1 r dw'w'{w-w')^,^u{w'euw'e2y-A-w') . (2.2) 
dOk w Jo 

Let now fj, be the magnification along the light ray, i.e., the solid-angle distortion of 
a light bundle. Then, {/j.) = 1 when averaged over the whole sky. The magnification // is 
determined by the determinant of 7^^, 

det{5jk+ljk) 

where it was assumed in the last step that |7jfc| <^ 1. The deviation of fi from unity, 
5 11 = IX — (ii) = fi — 1, is then determined by the trace of the shear tensor, or 

bii{d,w) = -bjklik{BM = — \ dw'w'{w-w')^ ' , ' -, 2.4 

w Jq (1 — w')^ 

where we have used Poisson's equation (All). Eq.(2.4) relates the magnification of the 
light ray to the density contrast along the light ray. 
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2.2 The angular QSO-galaxy correlation function 

Assume now that the observer sees a QSO population above flux threshold S charac- 

— * 

terized by its number density on the sphere nQ{S; 9), and similarly a galaxy population 

— * 

characterized by nG(^)- Their angular cross-correlation function is 

— * — * 

where the average extends over 9, and the result does not depend on the direction of 

due to the statistical isotropy of the QSO and galaxy distributions on the sky. 

The QSO population seen by the observer depends on the unlensed QSO number 

counts and, since the QSO sample is flux limited, also on the magnification of the QSOs 

via the magnification bias. Let the unlensed cumulative number density be 

n'^iS') (X (Sr , (2.6) 

where S' is the QSO flux in absence of gravitational magniflcation. The power-law 
index a is assumed to be independent of the QSO redshift, which appears to be an 
uncritical approximation. Due to gravitational lensing, the flux is magnified by a factor 
/J,. Therefore, S = [iS' , and the observed number counts are related to the unlensed 
counts by 

""f^) = I "^ if) ■ (^■^' 

The factor 1/fi out front accounts for the solid-angle distortion which causes the magni- 
fication; the solid angle of a magnified patch of the sky is reduced. With the power-law 
assumption (2.6), we then have 

nQ(5)ocAi-«-^5V (2.8) 

If a = — 1, the unlensed number counts are unchanged. Since the average number density 
of QSOs equals the unlensed one because {/j) = 1, we obtain 

'^q(^;^) -(^q('5)) _ ..-a-i 



{nqiS)) 



//-"-^ - 1 . (2.9) 



With 11 = 1 + S/d and 5/i ^ 1 as expected from weak lensing (see, e.g., Jaroszyhski et al. 
1990, Bartelmann & Schneider 1991), 

^-"-^ -l^{-a-l)5n. (2.10) 

We now assume that the galaxy number density is related to the density contrast 5 
by the bias factor b, 

^^^^^ = b5. (2.11) 

(«g) 

The efi'ect of a possible scale dependence of the bias factor b is addressed in the next 

section. Then, the correlation function Cqg{4>) of Eq.(2.5) reduces to 

eQG(<^) = {-a - 1)6 (^di,{9)d{9 + $)) = {-a - 1)6 (,,s{<l>) , (2.12) 

where ^^^ is the angular cross-correlation function between magnification and density 
contrast. 
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2.3 Evaluation of ^ns{4>) 

2.3.1 QSO weight function. We assumed that the QSO number counts have a flux dis- 
tribution with slope a independent of their distance. Their distance distribution be 
characterized by a normalized function Wq{w). Then, the magnification S/i of Eq.(2.4) 
has to be weighted with Wq, or 



(2.13) 



where 



5fi{9)= I dwWQ{w)Sfi{9,w) 
Jo 

Changing the order of integrations, we find 

Sm = 12 /' dw' GqK)^7?^^^-^ , (2.14) 

Jo {i-w'y 

/I 
dw' {l - ^^ Wq{w') . (2.15) 

If all observed QSOs were located at the same distance wq, we would have 

Wq{w) = 5{w - wo) , (2.16) 

which would result in 

Gq{w) = R{wo -w)w(l-—] , (2.17) 

where H(a;) is Heaviside's step function. 

Instead, we will adopt in the following the QSO weight function 

Wq{w) = Cw-^ , to €[100,1], (2.18) 

where Wq is the minimum distance of QSOs in the sample, which is related to the lower 

QSO cutoff redshift 

1 





^^ (1- 


way 












Wq to be normalized, we require 
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(2.19) 



(2.20) 



(2.21) 



where 

w = inax{w,wo) . (2.22) 
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Fig. 1. The weight functions Wq{w) for 
/3 = 1 and o; G {1, . . . , 4} (top left frame), 
for a = 1 and /3 £ {1,...,4} (top right 
frame), and Gq{w) for 7 £ {0, . . . , 3} in 
the bottom frame, zq = 0.2 and zq = 1.0 
are kept fixed 



2.3.2 Galaxy weight function. The galaxy distribution along the light ray must also be 
weighted by some galaxy weight function Wq{w) which is meant to mimic the distance 
distribution of galaxies in a magnitude-limited galaxy sample. We adopt the weight 
function given by Kaiser (1992), 



Wg{w) 






exp 



w 



(2.23) 



which is normalized to unity. Kaiser (1992) chooses a = 1 and /5 = 4. The average 
galaxy distance is related to w^ by 



wg = w^ 



r[{2 + a)/P] 

r[(i + «)//3] ' 



(2.24) 



and the average galaxy redshift zq is determined by wq analogous to Eq.(2.19). With 
Eq.(2.23), we obtain the projected density contrast 



m 



dw5{w, 1 — w)Wg{w) 



(2.25) 



Figure 1 shows the weight functions Wg{w) for /5 = 1 and a G {1,...,4} in the 
top left frame, for a = 1 and /5 G {1,...,4} in the top right frame, and Gq(w) for 
7 G {0, . . . , 3} in the bottom frame, zq = 0.2 and zq = 1.0 are kept fixed. 

2.3.3 Power spectra. Kaiser (1992) has shown that it is most convenient to evaluate 
correlation functions of the type considered here in Fourier space, making use of the 
fact that the auto-correlation function of some quantity p is the Fourier transform of the 
power spectrum of this quantity. 



/(Pk 
-——^Pp{K)exp{-iK(j)) 



(2.26) 



Both quantities of interest here, S/j, and 5, invoke weighted integrations over the 
density contrast along some light ray. The problem therefore arises of how the spectrum 
Pp{K.) of a quantity p which is a weighted projection of a three-dimensional quantity / 
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is related to the spectrum Pf{k) of /. The answer is given in the appendix of Kaiser 
(1992), Eq.(A2). Given a statistically homogeneous and isotropic random field f{x) and 
a weighted projection 

p{e)= I dw q{w)f{w0i,we2,w) (2.27) 

of /, then the spectrum of p is 

Pp{K)= f dw^-^Pf(-,l-w) . (2.28) 

Jo w^ \w J 

This Fourier-space analogue of Limber's equation (e.g., Peebles 1993, Eq.(7.47)) arises 
if it is assumed that the ^L'-integration can be decomposed into a sum over shells with 
width Aw <^ 1, whose contributions are mutually statistically independent such that 
Pf{k) ~ for A; ^ (l/Aw). For the purpose discussed here, this assumption is safely 
fulfilled. A typical scale in the power spectrum of the density contrast is ~ 50 Mpc/h, 
which is much smaller than twice the Hubble length which corresponds to w = 1. 

We are here after the cross-correlation function between magnification and density 
contrast ^^5. Both the magnification and the density contrast have the form of Eq.(2.27), 
where 12Gq{w)/{1 — w)"^ and Wg{w) adopt the role of q{w), respectively. In both cases, 
f{w) is the density contrast S{w). Therefore, we can insert Eqs.(2.14) and (2.25) into 
(2.28) and obtain 

P,, W = 12 /'d»52WI!°Mp, («,!_„). (2.29) 

Jo to^(l — to)^ \w ) 

The angular cross-correlation function ^^5 now follows from Eqs.(2.26) and (2.29); 

(2.30) 

— * 

Where we have introduced the angle ^ between k and ^. The rightmost integral results 
in the zeroth-order Bessel function of the first kind, 



I 



'^'^ dd 



27r 



— exp(— i/t^cosf^) = Jo(/«^) . (2.31) 

27r 

Further, we substitute k = (k/w) and obtain 

. .,^ f\ Gq(w)Wg(w) f^kdk^,^ ^. . , ,^ 

^^s{(p) = 12 / dw 7^ ^ ":; ' / -^Psik, 1 - w)M'^H) ■ 2.32 

Jo (1 - wY Jo 27r 

If we assume linear density evolution, we have 

Ps{k, l-w) = {l- w^Psik, 0) = (1 - w)^P^''\k) , (2.33) 

and thus 

e^,(<^) = 12 / ^—Pf\k) / dw{l - wfGci{w)WG{w)h{wk<l^) . (2.34) 

Jo 271 Jo 
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We will use for CDM the linearly evolved perturbation spectrum as given by Bardeen et 
al. (1986). For HDM, we will instead use the spectrum as derived from the Zel'dovich 
approximation by Schneider & Bartelmann (1994). Since the Zel'dovich approximation 
affects the spectrum only for small scales (large k), this spectrum will also be scaled 
in proportion to a^ = (1 — to)^ as in the linear case, but will be modified for large k 
dependent on 1 — w. 

To acquire some insight into the behaviour of ^^^(0) for CDM and HDM spectra, 
we insert the model spectra 

^i.HDM(^) = AuDukexpi-Xk) ; 
(0) _ 1 (2-35) 

into Eq.(2.34). For the exponentially decaying spectrum, the ^-integration is readily 
performed using Eq.(6.621.4) of Gradshteyn & Ryzhik (1980); it yields 

C^F''{4>) = ^HDM /' dwil - wrGc^iw)WGiw) ^f ~ ^if . (2.36) 

For the algebraically decaying model spectrum, the ^-integration leads to hypergeometric 
functions 1F2; the result is therefore of little practical use. However, series expansions of 
the correlation functions after evaluating the integral in Eq.(2.36) for both model spectra 
for small ^ yield 



i';:T{<^) = i?o - i^K/' + B^<\>^ + o{<\?) 



(2.37) 



where the constants Ci and Dj are 



and 



1 /-i 
Co = ^HDM— TT / dw K{w) , 

3 /"^ 

C2 = ^HDM— r^ / dw w^Kiw) , 

1 (^ 

Do = AcBM^r / dwK{w), (2.38) 

olio Jo 
1 [^ 

Di = AcDMTr / dwwK{w), 
2v5" Jo 

D2 = AcDM^ dww^K{w) , 
^^ Jo 

K{w) = {l-wfGci{w)WG{w). (2.39) 



Eq.(2.37) shows that, to second order in 0, the QSO-galaxy angular correlation function 
decreases like a parabola for HDM and approximately linearly for CDM. Figs. (2a, 2b) 



Numerical simulations 



0.001 



0.0001 

it 

0.001 



0.0001 




0.4 



4 2 

[degrees] 



Fig. 2a. The angular cross-correlation func- 
tion ^ij,s{(f>) for a CDM spectrum. The four 
frames of the figure show results for differ- 
ent choices of zq as indicated in the frames, 
and the line types distinguish between dif- 
ferent choices for zq as indicated in the top 
line 
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Fig. 2b. Similar to Fig.2a, but for an HDM 
spectrum 
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Fig. 3. Comparison between the angular 
cross-correlation function ^iis{4>) for HDM 
(left frame) and CDM (right frame), zq was 
set to 1.5, and the line types distinguish be- 
tween different choices for zq, as indicated 
at the top of the figure. The differences in 
shape and amplitude for the two density- 
perturbation spectra are clearly seen 



show ^jj,s{(l>) for CDM and HDM for various choices of zq (indicated by the hne type) 
and Zq (indicated in the four frames of these figures). 

Figure 3 directly compares the frames for zq = 1.5 from Figs. (2a, b) to emphasize 
the difference in amphtude and shape of ^^^(0) for HDM (left frame) and CDM (right 
frame) . 
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Figs. (2a, 2b, 3) show the quahtative behaviour of Cfis{(t>) derived in Eq.(2.37): The 
HDM correlation function starts with horizontal tangent at ^ = and declines slowly, 
while the CDM correlation function starts with a rather steep decline which flattens for 
larger 0. 



3 Numerical simulations 

This section describes numerical large-scale structure simulations performed to test 
whether the analytical angular cross-correlation function between QSOs and galaxies 
from gravitational lensing is reproduced numerically, to see how many QSOs one would 
need, and how deep galaxy observations should go in order to reliably estimate Cqg{4>)- 

To simulate the large-scale structure and the galaxy distribution associated with it, 
I use the technique described in detail by Bartelmann & Schneider (1992). Briefly, the 
large-scale structure is modelled using the adhesion approximation (Gurbatov, Saichev & 
Shandarin 1989, Weinberg & Gunn 1990). This proves useful because this approximation 
renders a very fast numerical code, which allows to simulate statistically independent 
dark-matter distributions on a large number of lens planes. The code provides a discrete 
mapping from Lagrangian to Eulerian space. To flnd the positions of 'galaxies', the 
density contrast along each trajectory starting from Lagrangian space is monitored in 
time. When it rises above a specifled threshold at some redshift for a given trajectory, 
this trajectory is assumed to carry a 'galaxy' later on. 

Speciflcally, the simulation starts from CDM initial conditions normalized to erg = 1 
in a flat universe (i7o = !> ^o = 0) with h = 1. The simulation boxes have a comoving 
size of (100 h~^ Mpc)^. Their redshifts are chosen such that the boxes are adjacent. The 
matter distribution in the boxes is projected onto their midplanes, yielding the surface 
mass density and thus the convergence required for lensing. 

The 'galaxies' are assigned luminosities drawn from the Schechter luminosity func- 
tion (Schechter 1976) with the parameters determined by Efstathiou, Ellis & Peterson 
(1988), 

u = -1.07 , M^ = -19.68 , if* = 1.56 x 10"^ Mpc"^ , (3.1) 

assuming h = 1. The Schechter function is cut off for luminosities below one tenth of L* , 
and their number density is chosen such as to agree with the normalization constant (f* . 
Fig. 4 shows the galaxy-galaxy autocorrelation function of the resulting galaxy distribu- 
tion. The figure shows that the autocorrelation function of the numerically positioned 
galaxies well reproduces the observed function. The magnitude-redshift relation for the 
model galaxies does not include /^-corrections or luminosity evolution. For further de- 
tails on the numerical large-scale structure model please refer to Bartelmann & Schneider 
(1992). 

Light rays are propagated through the simulated stack of lens planes using the 
multiple lens-plane approximation (see, e.g., Schneider et al. 1992, Chap. 9). Then, a 
synthetic QSO catalog is created assuming a random spatial distribution above redshift 
zq = 1, with redshift distribution such that the fraction of QSOs observed within dw 
of w is independent of w (corresponding to 7 = in Eq.(2.18)). Intrinsic fluxes S' for 
these QSOs are drawn from a cumulative luminosity function of the form (2.6), with 
a = —5. QSOs are included as visible in the catalog if iJ,S' exceeds a specifled flux 
threshold, where fi is the magniflcation factor along the line-of-sight to the QSO, and 
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Fig. 4. Galaxy-galaxy autocorrelation func- 
tion of the galaxy distribution created by 
the numerical simulation described here. 
The dotted line is the curve (r/ro)~ with 
ro = 7 Mpc/h 



Fig. 5. Angular cross-correlation function 
between QSOs and galaxies derived from 
the numerical large-scale structure model 
described here, divided by —a — 1=4. The 
curve has bootstrap error bars ranging from 
the 10- to the 90-percentile values found for 
each angular separation. The smooth curve 
shows the analytic result for Cqg('^)- Note 
that h = 1 was chosen here 



the solid-angle distortion due to lensing is taken into account by weighing the local QSO 
number density with /x~^. 

Prom this numerically created QSO catalog, 50 subsamples of 100 QSOs each are 
randomly chosen. From the distribution of galaxies with magnitudes < 21 in circular 
fields of 28' radius, the angular cross-correlation function between QSOs and galaxies is 
determined for each of the 50 subsamples. The average number density of galaxies re- 
quired for the normalization of the cross-correlation function is found by placing identical 
fields randomly on the simulated 'sky'. Fig. 5 shows the resulting angular cross-correlation 
function, divided by —a — 1 = 4. 

To compare the numerically determined cross-correlation function with the analytic 
result of Eq.(2.34), we have to know the bias factor of the simulated galaxy distribution. 
Since the galaxy-formation scheme employed here is not simple Eulerian biasing, the 
bias factor is not an input parameter of the simulation. Therefore, I determine b{r) as 
the ratio of the rms fluctuations in galaxy number density and matter density (see also 
Eq.(3.4) below), 






b{r) 



rms,r 



Sp_ 

(p) 



(3.2) 



this quantity b{r) is termed 'integral spatial bias' by Cen & Ostriker (1992) and Ostriker 



rz 



Losmological parameters from angular correlations between ybUs and galaxies 




Fig. 6. Integral spatial bias factor b{r) of 
the numerically simulated galaxy distribu- 
tion 



(1993). The so-determined b{r) is displayed in Fig. 6. 

The figure shows that the integral spatial bias factor b{r) exhibits the same qual- 
itative behaviour as derived by Cen & Ostriker (1992), in that it approaches ~ 2 at 
r ~ 1 M-pc/h and falls off towards unity for larger r. To incorporate a scale-dependent 
bias factor into the computation of the analytic angular cross-correlation function ^^^(0) 
of Eq.(2.34), we employ the differential spatial bias (Cen & Ostriker 1992) defined by 



b{k) 



Psik) ' 



(3.3) 



where -Pg(^) is the spectrum of the number-density fiuctuations of the galaxy distribu- 
tion. This quantity can be related to the integral spatial bias factor noting that the rms 
fiuctuations in both the galaxy number density and the matter density are given by 
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where W{kr) is the top-hat window function 

W{kr) 
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(3.5) 



with the spherical Bessel function of first kind ji. The symbol ax{r) abbreviates the rms 
fiuctuations of the quantity x (either the galaxy number density or the matter density) 
on the spatial scale r. For an approximation, note that the window function W{kr) can 
reasonably be approximated by the Heaviside step function H(l — kr). Then, making use 
of isotropy, Eq.(3.4) can be written 
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Upon differentiation with respect to r, we find from Eq.(3.6) 
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Fig. 7. Numerically determined QSO-gala- 
xy cross-correlation function ^qg(<^)i di- 
vided by — a — 1 = 4, for different depth of 
the synthetic galaxy 'observations' and con- 
stant QSO sample size A^qso = 100. The 
curves in the four frames were obtained for 
mmax & {19,20,21,22} as indicated 



where kr = (I/'')- Inserting this into Eq.(3.3) and using Eq.(3.2) yields 



b{kr) ~6(r). 



4TT^r^aj{r) d\nb{r) 



Psikr] 



dlnr 



(3.8) 



Since b{r) decreases with r, b{kr) > b{r). Having found b{k), we account for the scale- 
dependent bias factor by introducing the factor b{k) into Eq.(2.34), replacing PJ ' (k) 

by b{k)P^ (k)- Strictly speaking, b{k) can also depend on the conformal time (1 — w). 
However, since only galaxies at rather low redshifts contribute to the cross correlation 
between QSOs and galaxies, this possible time dependence is omitted. Finally, the dashed 
curve in Fig. 5 was determined this way, choosing 7 = in the QSO weight function 
Gq{w) of Eq.(2.18) to mimic the distance distribution of the synthetic QSOs in the 
numerical simulation, further choosing zq = 1, and a = 2, (3 = 1, w^ =0.3. The latter 
three parameters render a good fit to the galaxy distribution of the numerical simulation. 
As seen in Fig. 5, the analytic cross-correlation function fits the numerically determined 
function very well. Without accounting for the spatial dependence of the bias factor, the 
fit is still acceptable for 6 ~ 1.5 ... 2, but then the analytic cross-correlation function is 
slightly too shallow. 

To illustrate the effects of reducing the depth of the galaxy observations and of 
reducing the size of the QSO sample, I show in Fig. 7 the numerically determined QSO- 
galaxy cross-correlation function for mmax ^ {19, 20, 21, 22}, divided by —a — 1 = 4, while 
the QSO sample size A^qso = 100 is kept fixed, and analogously in Fig. 8 ^qg(0)/(— o — 1) 
for different synthtetic QSO sample sizes TVqso ^ {25, 50, 75, 100} and mmax = 22. 

From the close correspondence between the analytic and the numerical cross- 
correlation functions in Fig. 5, and from the changes in the cross-correlation function 
with mmax and iVqso shown in Figs. 7 and 8, I conclude that (1) the cross-correlation 
between background QSOs and foreground galaxies can indeed used as a tool to ob- 
tain direct information on the dark-matter density spectrum and, in particular, the bias 
factor, and that (2) this cross-correlation function can probably be measured observing 
galaxies down to ~ 2P* . . . 22"^^ magnitude in fields of about 20' . . . 25' radius around 
50 . . . 100 QSOs. Reducing the QSO number to ^ 50, or the limiting galaxy magnitude 
to mniax ^21, increases the error bars such that the signal can no longer be significantly 
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Fig. 8. Numerically determined QSO-gala- 
xy cross-correlation function ^qg(<^)! di- 
vided by —a — 1 = 4, for different QSO 
sample sizes and constant limiting galaxy 
magnitude rrimax = 22. The curves in 
the four frames were obtained for Afqso £ 
{25, 50, 75, 100} as indicated 



distinguished from zero. If the fields are reduced to much less than ~ 20' radius, the 
shape of the correlation function can hardly be determined, and much larger fields are 
hard to obtain in practice. 

The size of the error bars in Figs. 5, 7, and 8 are dominated by fluctuations in the 
galaxy number counts around the QSOs. It can be reduced by binning the galaxy counts 
in larger bins than used here. 

Of course, the result also relies on the QSO luminosity function being sufficiently 
steep, since, as shown in Eq.(2.12), the QSO-galaxy cross-correlation function is propor- 
tional to (a — 1). For a = —1, the result would be zero. However, I show in the next 
section that for the 1-Jansky sample of radio-loud QSOs, the double-waveband magni- 
fication bias renders effective values for a which are compatible with the value a = — 5 
chosen here. 



4 Double- waveband magnification bias 



It was first argued by Borgeest et al. (1991) that the magnification bias is more efficient 
in a sample of sources which is fiux-limited in two or more wavebands rather than in only 
one, provided the fluxes in these wavebands are mutually uncorrelated. The reason is easy 
to see. Consider a sample of sources which is flux limited in n different wavebands. The 
unlensed counts of these sources depend on n fluxes S'^. Assuming that the magniflcation 
factor in all n wavebands is the same, Eq.(2.7) reads 
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(4.1) 



The latter assumption that the magniflcation factor is the same in all n wavebands is not 
trivial despite the achromaticity of lensing because the radiation in different wavebands 
can be emitted from regions of different size of the sources. Microlensing, e.g., can 
only affect small parts of the source like the continuum-emitting region, while the radio 
emission usually comes from much larger regions whose flux remains unaffected. 

The double-logarithmic slope of the number counts effective for the magniflcation 
bias is given by 
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Fig. 9. Effective cumulative slope of 1-Jans- 
ky sources with redshift 2; > 1 as a function 
of the limiting optical magnitude rrimax- 
For optically bright subsamples, the effec- 
tive number-count slope decreases to Oeff = 
-4... -5 
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where Eq.(4.1) was used. As argued above, the magnification factors induced by the 
lensing effect of large-scale structures is close to unity, /j, = 1 + S/j, with 5// <C 1. There- 
fore, only the local behaviour of the unlensed number counts is relevant for the present 
consideration. If the number counts can locally be approximated by power laws with 
index q, 

d\nn'(S'...,S') 

and if the unlensed fluxes S'^ are mutually independent, it follows from Eq.(4.2) that the 
effective double-logarithmic slope is just the sum of the slopes in the n wavebands. 
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(4.4) 



In any case, Eq.(4.2) shows that imposing flux thresholds in two or more wavebands 
effectively steepens the number-count distribution. 

This steepening of the effective number counts is demonstrated in Fig. 9 using the 
optically identifled subsample of the 1-Jansky sample of radio-loud QSOs with redshift 
z > 1. These are by deflnition flux- limited to Sy-ad > 1 Jy at 5 GHz, and since they 
are optically identified, one can additionally impose an optical flux limit. For the figure, 
optical magnitudes m of the 1-Jansky sources are converted to fluxes relative to some 
arbitrary reference magnitude m^a,x, 

Sopt = loO-4("^-"^--) . (4.4) 

For S > 1, the cumulative slope a^fi of the counts Popt,rad('S') of 1-Jansky sources with 

min(5opt,5rad)>5 (4.5) 

is determined and plotted as a function of mmax- The number-count slope without optical 
flux limit is a^ad = —2.3. 

The flgure shows that subsamples of the 1-Jansky sources can have a steep effective 
number-count slope when they are flux-limited in the optical wavelength in addition to 
the intrinsic flux limitation in the radio waveband at 5 GHz. 
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5 Summary and discussion 

It was shown in previous papers that there is evidence for correlations on angular scales 
of several 10 arc minutes between background QSOs and foreground galaxies on a high 
significance level (Fugmann 1992, Bartelmann & Schneider 1993b, 1994a, Bartelmann 
et al. 1994b). As argued there and in Bartelmann & Schneider (1993a), gravitational 
lensing by large-scale structures provides a viable explanation for these correlations. If 
this interpretation is valid, it proves the biasing hypothesis of galaxy formation and 
provides a method to trace extended dark structures by their gravitational lens effect, 
which would otherwise be hard to detect. 

This paper develops this idea further. I have computed the angular two-point cross- 
correlation function ^qq between background sources and foreground galaxies due to the 
gravitational lensing effect of large-scale structures. This correlation function is sensitive 
to the power spectrum of dark-matter fluctuations and therefore potentially allows to 
discriminate between different cosmogonic models. As shown in Eqs.(2.12), (2.37), and 
(2.38), ^QG is proportional to the spectral amplitude A, to the bias factor b of galaxy 
formation, and to (—a — 1), where a is the double-logarithmic slope of the cumulative 
QSO number counts. The spectral amplitude A is known from the analysis of the CMB 
fluctuations, and the slope a can independently be determined from the background QSO 
sample employed. Thus, ^qg also allows to measure the bias factor b. 

Numerical simulations suggest that ^qg can be measured using a sample of 50 ... 100 
QSOs with redshift z > 1 ii galaxy observations down to a limiting magnitude of 21 ... 22 
are available in flelds of about 25' radius around these QSOs. The slope of the QSO 
number counts can effectively be increased by imposing flux limits in two or more wave- 
bands upon the QSO sample. This multiple-waveband magniflcation bias (Borgeest et 
al. 1991) can increase the effective number-count slope to a^a ^ — 5 in the 1-Jansky 
sample of high-redshift, radio-loud QSOs. 

There are several problems which may render the determination of ^qg difficult in 
practice. First, ^qg has to be properly normalized. To do so, the average galaxy number 
density has to be determined, a measurement prone to observational errors. However, 
since ^qg is an angular correlation function, it is much easier to measure than the spatial 
galaxy autocorrelation function, because errors in the distance determination of galaxies 
do not aff'ect the result. Even if the amplitude of ^qg is not determined correctly, its 
shape can still yield information about the dark-matter perturbation spectrum. Sec- 
ond, galaxy counts and optical QSO fluxes are influenced by dust absorption. Although 
the 1-Jansky QSO sample is radio selected, this effect could become important when 
one imposes an additional optical flux limit to take advantage of the multiple-waveband 
magniflcation bias. More severe, however, is the influence of spatially varying dust ab- 
sorption or reddening on the galaxy number counts. This could affect the normalization 
and the shape of ^qg- In order to minimize such effects, the galaxy counts required for 
determining ^qg should be done in red light or near infrared. Third, the correlation scale 
investigated here requires large CCD frames or CCD mosaics. While sufficiently large 
CCDs will probably become available in the near future, CCD mosaics pose the problem 
of adapting galaxy counts to a common level in all their parts. 

Despite these observational problems, the measurement proposed here provides a 
tool to determine the shape of the power spectrum of density fluctuations independent 
of other attempts, and to quantify the biasing hypothesis of galaxy formation; the latter 
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being hard to fulfil in a different way. 

Acknowledgements . I wish to thank Nick Kaiser and Peter Schneider for fruitful discus- 
sions and valuable comments. 

Appendix A. The deflection angle of a weakly perturbed light 
ray 

We restrict ourselves to a model universe which is on average of the Einstein-de Sitter 
type (i7o = 1, ^0 = 0, p = 0) and whose density is weakly perturbed by the relative 
density contrast S with (S) = when averaged over sufficiently large volumes. In close 
analogy to Kaiser (1992), we write the space-time metric in the form 

ds"^ = a^{r}) [dTf - Sijdw'-dw^] , (^1) 

where rj is the conformal time which is defined by 

c dt c da , , , 

dr] = —— = —^ (A2) 

a{t) aa 

and the condition that r; = when i = 0. Friedmann's equation reads, for our choice of 
the cosmological parameters, 

a = ^ , (^3) 



a 

and since a{z) = {1 + z)~^ when a is normalized to unity at present, ao = a{to) = 1, we 
find from Eqs.(A2) and (A3) 

Ho y/l + Z 

Thus, the comoving time has the dimension of a length. If we express in the following 
all lengths in units of twice the Hubble length, we can write 

viz) = -1^ = v^ • (^5) 

The comoving distance w along a radial light ray starting at the observer position 
follows from ds = 0, or, from Eq.(Al), 

drj = —dw , (^6) 

and therefore the comoving distance equals the conformal lookback time, 

l-r^{z) = w{z) = l--^. {A7) 

The density contrast S will give rise to gravitational lens effects on light rays. It 
satisfies Poisson's equation, which reads, in comoving coordinates, 

A$ = Ai^Gpa^S ; (^8) 
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see, e.g., Peebles [1993, Eq.(5.107)]. p is the mean density of the background universe, 
hence 



and thus Eq.(A8) transforms to 



otj2 

2a ^ ' 



Measuring lengths in units of twice the Hubble length and the potential in units of c^, 
we can write instead of Eq.(AlO) 

A$iw,i]) = -5{w,T]) = —5{w,T]) , (^11) 

a r]^ 

where we have inserted the arguments for clarity: w is the comoving position considered, 
and 77 is the conformal time. 

Consider now a radial light ray starting at the observer position into direction 
(^1,^2); then, w = {w9i,w92,'w) and 77 = 1 — w. This light ray will be deflected by 
the gravitational-lens effect of the density perturbations according to the equation of 
motion 

rj{w,r]) = -2$j{w,l-w), {A12) 

where r is the comoving deviation vector of the light ray relative to the unperturbed ray, 
and the comma preceding the index j denotes the partial derivative with respect to rj. 
The total deviation from the unperturbed ray is then given by 

pw pw 

rj{ei,e2,w) = -2 / dw' / dw"^j{w",l-w") . {A13) 

Jo Jo 

This inner integral is easily evaluated along the unperturbed light ray; it yields 

pw 

rj{ei,92,w) = -2 dw'{w-w')$.j{w'di,w'e2,w';l-w') . (AU) 

Jo 

The deflection angle of the light ray propagating into the direction 9 = {61,62) is given 
by waj = rj, or 

aj{e,w) = / dw'{w-w')^^j{w'ei,w'e2,w';l-w') . {A15) 

^ Jo 

This is Eq.(2.1) for the deflection angle of a weakly perturbed light ray propagating 
through an Einstein-de Sitter background universe. 
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